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Friction between ring-polymer brushes at melt densities sliding past each other are studied using 
extensive course-grained molecular dynamics simulations and scaling arguments, and the results 
are compared to the friction between linear-polymer brushes. We show that for a velocity range 
spanning over three decades, the frictional forces measured for ring-polymer brushes are half the 
corresponding friction in case of linear brushes. In the linear-force regime, the weak inter-digitation 
of two ring brushes compared to linear brushes also leads to a lower number of binary collisions 
between the monomers of opposing brushes. At high velocities, where the thickness of the inter- 
digitation layer between two opposing brushes is on the order monomer size regardless of brush 
topology, stretched segments of ring polymers take a double-stranded conformation. As a result, 
monomers of the double-stranded segments collide less with the monomers of the opposing ring 
brush even though a similar number of monomers occupies the inter-digitation layer for ring and 
linear-brush bilayers. The numerical data obtained from our simulations is consistent with the 
proposed scaling analysis. Conformation-dependent frictional reduction observed in ring brushes 
can have important consequences in non-equilibrium bulk systems. 


INTRODUCTION 

If polymer chains are grafted by one of their ends to 
a planar or curved surface, above a certain critical graft- 
ing density p* « 1/B% [I] , where Ro is the characteris¬ 
tic equilibrium chain size, the chains are stretched away 
from the surface due to steric repulsion from surround¬ 
ing chains and form polymer brushes. A polymer brush 
is a soft polymeric material that can deform under var¬ 
ious external forces. The external force can be due to a 
fluid flowing over the brush, a flow due to the relative 
motion of a second brush, or alternatively, an external 
electrical field (if the polymers are charged). Once the 
external stimulus is completely removed, as the chains 
constituting the brush relax, the deformation of brush 
is reversed similar to the elastic deformation observed 
for solids. However, the tribological behaviour of poly¬ 
meric systems more resembles that of a fluid rather than 
a solid [2]. For instance, if two inter-digitated polymer 
brushes are slid past each other, frictional forces due to 
the relative motion of brushes vanish linearly as the rela¬ 
tive velocity of the brushes is decreased towards zero. In 
other words, the friction is viscous, for which no static 
friction occurs, unlike the solid-state friction, where finite 
forces are needed to initiate motion [3]. 

The polymer brushes has attracted increasing atten¬ 
tion due to their applications in nanotechnology and 
material sciences as bio-sensors, bio-fueling, stimuli- 
responsive surfaces ms, or for the stabilisation of col¬ 
loidal solutions mm- The most fascinating application 
of brush-like structures is facilitated by nature in main¬ 


taining the lubrication in tissues fTTHH]. For instance, 
in mammalian joints, where very low lubrication should 
be retained under pressures as high as 5 MPa, brush-like 
structures in combination with the synovial fluid provide 
lubrication in between the articular joints [THE]. The 
surfaces separating the articular cartilage and the syn¬ 
ovial fluid are thought to be covered by high molecular 
weight molecules such as lubricin min]. In turn, these 
long and charged glycoproteins may function as water- 
based bio-lubricants. The relationship between morphol¬ 
ogy of these long and charged macromolecules on the 
cartilage surfaces and their function in lubrication is still 
under investigation mmmm- Amphiphilic lubricin 
chains can adsorb on both hydrophilic and hydrophobic 
surfaces naEB, but their conformations in the adsorbed 
state strongly depend on the type of surface. These 
molecules can bind onto hydrophilic surfaces via their 
charged domains located in the central part of a molecule 
and form structures resembling a brush composed of lin¬ 
ear chains. Alternatively, they can also bind onto hy¬ 
drophobic surfaces with their terminal groups to form 
loop brush-like structures OEUI22]. Hence, based on 
experimental observations and variations in the molec¬ 
ular conformations in vitro , one may conclude that in 
reality the whole cartilage surface resembles a polymer 
brush composed of linear and loop-like chains. Although 
molecules mentioned so far are highly charged, and long- 
range interactions might be a dominant factor in the re¬ 
duced frictional forces, conformation of chains - whether 
they are linear or looped - can influence inter-digitation 
of molecules as well as their relaxation times, and hence, 
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friction and lubrication in tissues. 

Inter-digitation effects due to ring topology have also 
been reported in the context of genetic material: in a 
computational study, unconcatenated chains of ring poly¬ 
mers showed a weak trend towards inter-mixing with each 
other compared to mixing behaviour of linear chains in 
confined environments [23] . We should also underline 
that ring polymers in bulk and at melt densities exhibit 
different behaviour with respect to their linear counter¬ 
parts [241426] . The behaviour of ring polymers cannot be 
described by one single fractal dimension unlike linear 
chains in their melts: while size of a linear chain at melt 
state scales as Rol ~ N 1 / 2 , a ring polymer in the melt of 
rings is much more compact and the characteristic size 
scales Ror ~ TV 1 / 2 for short chains and Ror ~ TV 1 / 3 for 
sufficiently large polymerization degrees, N. In between 
these two limits, an intermediate regime Ror TV 2 / 5 
arises [24] [25]. Moreover, even in the regimes where the 
size of ring chain is non-Gaussian, interestingly higher 
moments of the end-to-end distance exhibit a Gaussian- 
like behaviour m- 

Existing computational [27ff33] and theoretical stud¬ 
ies [29] [341437] on linear-polymer brushes have shown that 
frictional forces acting brushes exhibit a cross-over from a 
linear to a non-linear regime upon increasing shear veloc¬ 
ity. The onset of non-linear friction typically appears at 
around shear rates where grafted chains begin to stretch. 
Ideally, one would expect a similar behaviour in the case 
of ring-polymer brushes. However, how the topology and 
peculiar scaling of ring chains alter the friction is an open 
question. To the best of our knowledge there are very 
few computational investigations of ring or loop poly¬ 
mer brushes out of equilibrium [38] [39] . In a previous 
work [39], tribolgy of loop brushes near the overlap con¬ 
centrations were studied, but the observed difference in 
frictional forces was negligible. However, at melt den¬ 
sities where the correlation length - distance between 
two chains - is on order of monomer size, the situation 
can be different as the number of collision between brush 
monomers can significantly increase in a denser system. 

Characterizing the systematic reduction in brush fric¬ 
tion due to the topology of the constituting chains can 
improve our understanding of how nature handles fric¬ 
tion and help the design and improvement of advanced 
biomimetic lubricants. Thus, spurred by the abundance 
of brush-like structures in nature and the interesting 
nature of ring polymers themselves in this paper we 
aim to study non-equilibrium behaviour of ring-polymer 
brushes at melt densities. In our extensive coarse-grained 
MD simulations, we used neutral (uncharged) polymer 
brushes. Given the fact that the system we would like to 
mimic is under high pressure (e.g., in joints), and ionic 
condensation (short Debye length in physiological con¬ 
ditions) can effectively neutralize the chains, we believe 
that this approximation is reasonable for the sake of min¬ 
imizing computational cost [28] since long-range electro¬ 


static interactions are known to be computationally ex¬ 
pensive, particularly for dense systems. Through a de¬ 
tailed analysis of simulation trajectories and by employ¬ 
ing scaling arguments, we demonstrate that the topology 
of chain in a melt brush is an important factor in the 
reduction of frictional forces. 

In our simulations, we found that for untangled 
brushes friction forces between two brushes made of lin¬ 
ear chains are always roughly a factor of two higher than 
those for ring chains. Although this difference is small, 
it persists for various grafting densities and chain sizes. 
This difference in the frictional forces of ring and linear 
brushes can be elucidated by the size of the overlap zone, 
which defines the amount of inter-digitation between two 
brushes. It turns out that the low tendency of ring chains 
to overlap with the opposing ring-brush leads to lower 
friction forces whereas linear chains can diffuse through 
the opposing brush more easily. Hence, brushes made of 
linear chains exhibit higher frictional forces. At very high 
velocities, on the other hand, segments of a ring brush 
take a double-stranded conformation. In turn, double¬ 
strands are less efficient in momentum transfer between 
the opposing brushes. The difference in friction forces 
between both systems is confirmed by scaling analysis. 

The paper is organized as follows: First we briefly de¬ 
scribe the simulation methodology. Next, the results ob¬ 
tained from non-equilibrium coarse-grained brush simu¬ 
lations will be discussed. By carefully analysing our data, 
we will relate the difference in friction forces to the in¬ 
trinsic properties of ring and linear (grafted) chains. The 
scaling arguments for brush systems will also be discussed 
to infer the difference in forces along with the simulations 
data. We conclude the paper by the summary of our find¬ 
ings and future prospects. 


SIMULATION DETAILS 


Simulations of polymer brush bilayers were performed 
using coarse-grained Kremer-Grest (KG) bead-spring 
model [31] [40]. Each individual chain of a polymer brush 
was composed of N (or 2 N) monomers (beads) con¬ 
nected by bonds. Two adjacent effective chain monomers 
are bonded via a “Finitely Extensible Nonlinear Elastic” 
(FENE) potential 


yF ENE( r ) = _l fcr 2i n x 



( 1 ) 


where bond stiffness k = 30e/cr 2 , distance between the 
beads r = |r|, and maximum bond length ro = 1.5 a m • 
The interaction strength e was measured in units of ther¬ 
mal energy ksT. The non-bonded interactions between 
all monomers were modeled by the truncated and shifted 
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Lennard-Jones (LJ) potential 

4e [(a/r) 12 — (a/r) 6 + 1/4] r < r c 
0 r > r Cy 

(2) 

where a was chosen as the unit of length and r c is cutoff. 
For polymer-polymer interaction we have used e = 1 and 
r c = 2 1 / 6 a. This choice of LJ potential cutoff in combi¬ 
nation with a monomer density of p m « 0.6 a -3 provides 
correct melt statistics mmm- 

Polymer brushes composed of chains with either linear 
or ring-like topology were studied. We performed simula¬ 
tions with linear-polymer brushes composed of N = 60, 
100 monomers, and ring-polymer brushes of N = 120, 
200 monomers per grafted chain. Chains were grafted 
on a square-lattice surface with dimensions 42cr x 36 a. 
The ring chains were grafted by one of their monomers 
on the surface. No equation of motion was solved for sur¬ 
face monomers (surface monomers were immobile during 
simulations). The surface and the brush monomers in¬ 
teract via Eq. ([!]), and Eq. (J2|. The grafting densities of 
linear chains are p R = 0.11 a -2 ,0.25a -2 , which respec¬ 
tively corresponds to inter-anchored monomer distances 
of 3 a, 2a at the surface. In the case of ring brushes, to 
obtain an equal monomeric density, p m ~ 0.6 a -3 , to 
the linear-brushes, the grafting densities were taken as 
p R = 0.5 p^ at the same inter-plate distance D. Note 
that D and p R ' R are connected via pm « N p R,R /D. In 
order to construct brush bilayer systems, first two non¬ 
interacting single brushes were generated as mirror im¬ 
ages of one another. While one of the brushes is fixed 
at z = 0, the other brush is brought into contact slowly 
at the desired inter-plate distance z = D to obtain the 
same monomer density for each p R,R and N. Finally, the 
systems are run at v = 0 for at least 10 7 MD steps to 
allow chains to relax. 

The molecular dynamics simulations were performed 
by solving the Langevin equation of motion, which de¬ 
scribes the Brownian motion of a set of interacting 
monomers, as 

mri = F EJ + F EENE - Tvi + Ff, i = 1,... , IV, (3) 

where = [xi,yi,Zi] is the position of i-th monomer. 
F/ J and Ff ENE in Eq. Q are respectively LJ and FENE 
forces exerted on the i-th monomer and given by deriva¬ 
tives of Eqs. 0 and 0 with respect to r^. The effect of 
the implicit solvent in Eq. (J3| is split into a slowly evolv¬ 
ing viscous force —Tr^ and a rapidly fluctuating stochas¬ 
tic force Ff . This random force Ff is related to the fric¬ 
tion coefficient T by the fluctuation-dissipation theorem 
(Ff(t)Ff(P)) = kBTTSijS(t — t'). The friction coefficient 
used in simulations was T = 0.5 mr" 1 , where m = 1 is 
the monomer mass and time was measured in units of 
r = /ma 2 /e. The integration step was taken to be 
At = 0.002 r. The velocity Verlet scheme was used for 



FIG. 1: Scheme of brush bilayer shear simulation. The 

force / acts in the direction opposite to velocity v. 

Chains at the top and bottom brush are rendered in 
different colors. The overlap zone (OZ) with thickness 
$l,r is indicated by a milky-white region. The dashed 
circle shows the segments of chains with gL,R monomers 
inside the penetration zone. Here N = 60, 
p R = 0.25 a -2 and D = 50 a. 

numerical integration of equations of motion Eq. 0- All 
simulations were performed in the N m VT- ensemble, i.e., 
constant volume V, total particle number N m and tem¬ 
perature T. The system temperature was set to the value 
T = 1.68 e/ks with fc# = 1 [41]. Periodic boundary con¬ 
ditions were introduced in the lateral directions, i.e. in x 
and y whereas in the ^-direction fixed boundary condi¬ 
tions were imposed. Simulations were carried out using 
molecular package LAMMPS m ■ Simulation snapshots 
were rendered using Visual Molecular Dynamics (VMD) 
[43]. 

The non-equilibrium shear simulations were performed 
as shown in Fig. [l] Both plates grafted by polymers 
were moved laterally in the opposite (±x-directions) at 
prescribed velocities in a range of v ~ 10 _3 a/r and 
v ~ 1 a/r. The inter-plane distance D was kept constant 
during shearing. For each plate velocity, all system were 
run for 10 7 MD steps until the steady state in friction 
force was reached, i.e. error bars in time-averaged quan¬ 
tities were independent of simulation time. Error bars 
were calculated via block-averaging. For the purpose of 
data analysis, additional simulations were run for 10 6 up 
to 10 8 MD simulation steps depending on the velocity 
and in order to obtain proper plate displacements. As 
a rule each plate was displaced by at least 5 times in 
the corresponding zbx-directions. To avoid any bias on 
friction forces or vertical chain diffusion while the sys¬ 
tem was sheared, the thermostat was applied only in the 
^-direction ED EH. 

In this paper, we report frictional forces that are the 
forces acting on the plates to keep them at a constant 
velocities, as well as averaged brush properties that are 
calculated from monomer trajectories by employing in- 
house analysis programs. Unless otherwise noted, all re- 
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suits presented in this paper were averaged over time. 


RESULTS AND DISCUSSION 

In the course of simulations, two compressed polymer 
brushes were moved in the opposite (=br) directions by 
driving both planes grafted by chains at prescribed ve¬ 
locities v as shown in Fig [l] Each plane experiences a 
friction force in the opposite direction to its motion due 
to relative motion of brushes. In Fig.[2j the friction force 
for linear brushes /l (filled symbols) and for ring brushes 
/r (open symbols) are shown as a function of velocity, 
u, for various brush systems with different polymeriza¬ 
tion degrees N and chain topologies. The overall veloc¬ 
ity dependence of friction forces is consistent with pre¬ 
viously reported results [28H3Q1 144] : The frictional force 
increases linearly up to a threshold velocity, which occurs 
at different velocities for different brush bilayers. Above 
a threshold velocity, a sublinear increase occurs for all 
brush systems considered here. We observe this linear- 
to-sublinear transition in most of our systems except the 
case with inter-plate distance D = 35 cr, where only the 
onset of the non-linear regime can be seen in Fig. [2ji. As 
we will discuss in more details in the following sections, 
this is due to the fact that the threshold velocity sepa¬ 
rating the linear and non-linear force regimes depends on 
segment size inside the overlap volume, where two oppos¬ 
ing brushes can co-exist. Hence, as the size of the aver¬ 
age segments increases inside the overlap volume, much 
slower velocities are required to observe the linear-force 
regime. The comparison of friction forces acting on lin¬ 
ear and ring brushes, which is the main motivation of 
this work, shows that friction forces for linear brushes 
are always higher than those acting on ring brushes, i.e. 
/l > /r- The numerical value of the ratio between these 
two forces fh/fR ~ 2 as can be seen in the insets of Fig. [2j 
The ratio holds for a broad range of shearing velocities 
(v ~ 10 _3 -10° cr/r) used in simulations (see the insets of 

Fig. [ 2 ). 

Average normal forces acting on the plates in the z- 
direction were also measured: We observed that normal 
pressures exhibit a weak dependence on velocity, i.e. a 
three order of magnitude increase in velocity leads to a 
roughly 10% increase in the normal forces for both linear 
and ring brushes (data not shown). At high velocities, 
a slight decrease in all measured pressures is observed. 
Similar to the frictional forces, the normal pressures mea¬ 
sured for linear brushes are factor of two higher than ring 
brushes, /p^ ~ 2. However, if total normal forces are 
rescaled by the number of grafted chains for each system, 
normal forces per chain are of almost equivalent in both 
systems. 

In what follows, we will discuss the observed difference 
in the frictional forces acting on linear and ring brushes. 
We will consider linear and non-linear force regimes sep- 



V [g/t] 


FIG. 2: Total friction forces / as a function of plate 
velocity v for linear and ring brushes. The grafting 
densities are a) p^ = 0.25cr -2 b) p^ = O.llcr -2 c) 

Pg = 0.25cr -2 d) p^ = O.llcr -2 with p^ = 0.5 p^ for all 
plots. 


arately and demonstrate that the ratio /l//# ~ 2 is 
related to the topology of chains and of their velocity- 
dependent conformations. 
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Linear regime 


In this subsection we will consider the regime where 
frictional forces increase linearly with the driving veloc¬ 
ity. If two dense polymer brushes are brought into con¬ 
tact at the inter-plate distance D > Rq, as shown in 
Fig [I] where Ro is equilibrium size of a free chain in bulk, 
chain segments near the free ends of grafted chains can 
interpenetrate through the opposing brush. Hence, the 
overlap zone (OZ) where monomers of opposing brushes 
can mutually interact can be defined (see a rectangular 
region in Fig. [I]). 

From simulation trajectories, the thickness of the OZ 
can be obtained using the cross product of monomer den¬ 
sity profiles of top p top (z) and bottom p hott (z) brushes. 
The cross product is non-zero only if monomers from 
both parts of brush coexist at the same z coordinate. 
The width of the OZ can be calculated from the cross- 
pro duct-weighted averages as 


S(v) = 2 




where we defined the cross product as uj(z,v) = 
Pto P (z,v)p bo tt(z,v). Hence, from Eq. Q, equilibrium 
value of OZ thickness reads S(v 0) = 5o- Note that in 
this paper we use subscript “zero” to indicate equilibrium 
(or linear-response) values of corresponding parameters 
(i.e. plate velocity v = 0) whereas subscripts L and R 
refer respectively to linear and ring brushes. 


The equilibrium values of the widths of the OZ for lin¬ 
ear and ring brushes Sol,or are shown in Fig.[3]as a func¬ 
tion of the number of monomers goL,OR per segment of a 
chain inside the OZ for linear and ring brushes. As a rule 
a segment is considered as being present inside the OZ if 
z coordinate of the TV-th monomer of the corresponding 
chain satisfies the following condition: D/2 — Sl,r < z < 
D/2 + Sl,r- As we perform our simulations nearly at 
melt density (p m ~ 0.6 cr -3 ) the random-walk statistics 
allows us to express square-root of the end-to-end dis¬ 
tance tol,or = (r’oLOi?) 1 ^ 2 °f a segment inside the OZ 

and composed of g 0 L,0R monomers as r 0LfiR ~ g^oR- 
The latter expression is valid for short segments. If the 
size of the segment determines the width of the OZ, then 
the relation vol,or ~ Sol,or should also hold. Indeed, 
this is what we observe in Fig. [3] and its inset for both 
linear and ring brushes: r 0 L,0R ~ S 0 l,or ~ 9ol,or- Inter- 
estingly, we found that at the same inter-plate distance 
D grafted chain in ring brush can occupy the overlap vol¬ 
ume with more monomers than a chain in linear brush 
(:qr > g L ). This finding demonstrate that inside the 
OZ segments of ring brushes are more compact as com¬ 
pared to segments of linear brushes. This fact is also 
supported by Fig. [3] since the pre-factor of the power- 



FIG. 3: The equilibrium values of the width of the 
overlap zone (OZ) versus number of monomers go per 
segment of a chain located inside the OZ for linear (full 
symbols) and ring (empty symbols) brushes. Dashed 
lines represent fitted power laws: So ~ 1.0 g^j 2 (for 
linear chains) and ~ 0.62hg 0 (f° r rings). The inset 
shows fluctuation of the end-to-end size of segments 
ro = (r^) 1 / 2 inside the OZ as a function of S q. Solid line 
stands for fitted scaling ro oc 5o- 


law fit is 1.6 times larger for linear brushes. The com¬ 
pactness of segments in ring brushes yields to a narrower 
OZ. As a remark, for very long segments in ring brushes 
{gR 1) the width of the OZ should in the asymptotic 
limit converge to ( r o r ) 1/2 ~ $ 0 R ~ gl R ■ In that case: 
Sol/Sor ~ fi'of 2 /% R - 

Since two opposing brushes can only interact inside the 
OZ, the observed frictional forces should be related to the 
width of OZ. The grafted chains are not static as they 
constantly diffuse in and out of the overlap volume which 
is given by A x Sl,r where A is the area of a grafting plate 
due to thermal fluctuations. When two plates are moved 
at v > 0, any chain segment that enters the OZ feels a 
flow induced by the relative motion of monomers mov¬ 
ing in the opposite direction. The force acting on each 
monomer inside the overlap volume can be most gener¬ 
ally expressed via Stokes drag f m ~ (/u, where ( is the 
monomeric friction coefficient. The total friction force 
acting on pm AS monomers inside the overlap volume can 
be expressed for both linear and ring brushes as follows 

/l,r ~ PmS l^rALIj^ r. ( 5 ) 

In Eq. ([5|, we defined Ql,r which quantifies the number 
of binary collisions between the monomers of two oppos¬ 
ing brushes. As we will see shortly, although there are 
A x S monomers inside the overlap volume, only a frac¬ 
tion of them participate in momentum exchange between 
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opposing brushes. Hence, Ql,r is an important quantity 
to distinguish the friction among different brush systems. 
In addition, the equilibrium value of the monomeric fric¬ 
tion coefficient (o is equal for monomers of both linear 
and ring brushes. This is due to the fact that on the time 
scales comparable with time between inter-bead collisions 
monomer cannot know whether it belongs to a linear or 
ring chain. 

Based on the linear-response theory, which states that 
any quantity takes its equilibrium value if the pertur¬ 
bation is small, we replace the quantities appearing in 
Eq. <§ with their equilibrium values as u —>• 0, i.e., 
Sr ~ Sol and Sr « Sor and similarly Qr « ^ 0 l and 
Qr ~ tioR- Thus, the frictional force ratio in the linear- 
force regime reads 

f_L ^ SpL^OL 

fR SoR^OR 

The validity of Eq. |6| can be verified from the simula¬ 
tion data by calculating the equilibrium values of Sl,r 
and &l,r for both linear or ring brushes. From now on, 
to keep the manuscript more compact we will only show 
data from four of our eight brush systems. All compar¬ 
isons and discussions are valid for the systems which are 
not shown here as well. We show our non-equilibrium 
simulation results for the thickness of OZ as a func¬ 
tion of plate velocity v for various linear and ring brushes 
in Fig.[4j Each frame of Fig. [4] compares the OZ width of 
ring-brushes (open symbols) to those obtained for linear 
brushes (filled symbols) at the same inter-plate distance 
D. Independently of the chain topology the width of the 
OZs exhibit a plateau for all brush systems as v 0. 
The plateau values of the OZs, which are Sl ~ Sol 
and Sr ~ Sor , persist up to roughly the plate veloci¬ 
ties where the linear-force regime observed in Fig. [2] also 
ends. Indeed, as we proposed for Eq. [6j the widths of 
OZs Sl,r are equal to their equilibrium values. This is 
due to the fact that within the linear-force regime, for 
which the frictional forces change linearly with velocity, 
values of Sol,or are nothing but the size fluctuations in 
the ^-direction of the chain segments inside the OZ. As 
also illustrated in Fig. [l] the size of a segment with ql,r 
monomers defines the width of the OZ. 

The number of monomers per segment 9l,r inside the 
OZ for linear and ring brushes was analyzed and is shown 
in Fig. [5] as a function of the plate velocity v using the 
same color code also as in Figs. [2] and [4j In the entire 
velocity range considered here the following ratio holds 
9r ~ 2 gL for the same value of D. At slow velocities 
(v 0) the ratio of monomers per segment gohjgoR is 
« 0.6. The ratio drops to Pl! 9r ~ 0.4 at higher velocities 
(v > 10 -1 cr/r). 

Similar trajectory analysis was also conducted to de¬ 
termine the number of binary collisions £2l,r inside the 
overlap volume: to count collisions, distances between 
the monomers from opposing brushes are calculated. If 



FIG. 4: The nonequilibrium thickness S of the overlap 
zone for linear and ring brushes plotted as a function of 
plate velocity v. The values of S were calculated using 
Eq. The insets show the ratios. The grafting 
densities are a) p^ = 0.25cr _2 and b) p^ = 0.11a -2 with 
p^ = 0.5Pg for all plots. 
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the distance is equal or smaller than the cutoff distance 
of LJ potential defined in Eq. ©> this specific pair is 
counted as a colliding pair. To obtain Ql,r and compare 
the number of collisions for ring and linear chains, the to¬ 
tal number of collisions is rescaled by the corresponding 
value of p m ASL,R- In Fig. [6] we show the binary-collision 
fractions Ql,r, as well as their ratios in the insets. In¬ 
terestingly, for linear brushes there are more collisions as 
compared to ring brushes, i.e. Qr > ^r for the entire 
velocity range although g r « 2 g^. An average ratio of 
Pl/pr ~ 0.5 indicate that even though a single chain in 
ring brushes can occupy the overlap volume with more 
monomers than a chain in linear brushes (gn > gL) the 
compactness of the segments in ring brushes yield to a 
narrower OZ. This is demonstrated in Fig. [3] where the 
pre-factor of the slope is bigger for linear brushes. The 
compactness of segments in ring brushes also results in 
a smaller number of collisions between the monomers of 
opposing brushes. Hence, /l//r ~ 2. In addition, as 
expected values of f \r and Qr are close to their values as 
v 0, which confirms that Qr ~ ^lo and CIr ~ Qo r. 

At the same plate separation D, the ratios of plateau 
values are Sqr/Sor ~ 1.3 db 0.05 from Fig. [4] and 
^ol/^or ~ 1.7 ± 0.1 from Fig. [6j Plugging the later 
expression into Eq. [6] gives the ratio of friction force for 
the linear-force regime Jr//r ~ 2.3, which is consistent 
with the ratios of forces obtained in the insets of Fig. [2|. 

As the frictional forces increase linearly (see Fig. [5]) in 
principle none of terms in Eq. © should have any veloc¬ 
ity dependence. Indeed, this is what one would expect 
at the linear-response level: at slow enough velocities 
(< 10“ 1 cr/r) the conformation of chain segments within 
the OZ are not affected by the velocity. As a result, the 
ratios of forces for linear and ring brushes are constant 
and given by Eq. <©• 

Our analyzes confirm that in ring-brush bilayers, the 
amount of inter-digitation of chains is weaker compared 
to that in linear-brush bilayers at slow driving velocities, 
where the friction increases linearly with velocity. 


Non-linear regime 

In this subsection we will focus on the regime where 
frictional forces demonstrate sublinear plate velocity v 
dependence. As can be seen in Fig. [2] nonlinear regime is 
observe for both linear and ring bilayers. The indication 
of the non-linear regime is the shrinkage of overlap vol¬ 
ume with increasing velocity v as shown in Fig. [4] Onset 
of the non-linear regime occurs at a threshold velocity 
v = v*. Indeed, v* has different values for different chain 
sizes as can be noticed in Figs. [2] and [4] Before we further 
continue our discussion on the non-linear regime, we will 
briefly discuss two scaling predictions for the threshold 
velocity, and refer them as v} and v}j. 

We will temporarily drop the indexes for ring and lin- 
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FIG. 6: The binary-collision factor inside the overlap 
zone between two brushes as a function of the plate 
velocity v. See the text for details. Insets show the ratio 
of linear-brush and ring-brush data. The grafting 
densities are a) p^ = 0.25 a~ 2 and b) p^ = O.llcr -2 with 
py = 0.5 p^ for all plots. 


ear brushes for simplicity and refer go = goR = 9or and 
So = Sol = Sor , etc. If a segment with go monomers en¬ 
ters the OZ, the force acting on this segment in the linear 
regime can be expressed as f c h ~ Co$o^- As we have dis¬ 
cussed in the previous section, below the threshold veloc¬ 
ity the chains are still Gaussian. Hence, as expected from 
a Gaussian chain, segment sizes in x, y and indirections 
are not coupled, i.e. fluctuations in size in all direction 
are independent of each other. However, as the velocity 
is increased, the frictional force per segment inside the 
OZ reach a threshold force f* h « ksT/b « (ogov} at v = 
vj [35] . The force on the segment given by f* h ~ ksT/b 
is high enough to align individual bonds in the direction 
of the relative motion. Thus, vertical segment fluctua¬ 
tions in the ^-direction, which define the width of OZ, 
will be affected by the shear in the acting in ^-direction. 
As a result the OZ decreases with respect to its equilib¬ 
rium value S < £o at v > v} ~ ksT/b(o9o ~ l/#o [35] . 
On the other hand, in a concentrated solution one may 
argue that the threshold velocity occurs at slower veloc¬ 
ities and is given by f* h ~ Co gv* n ~ ksT/So that leads 
to Vjj ~ ksT/bCog^/R [29] . where v is the scaling ex¬ 
ponent (v = 1/2 for an ideal chain and v = 0.588 for a 
swollen chain in 3D). In our simulations, hydrodynamic 
interactions are screened and we take v = 1/2 and ob- 




















tain Vjj ~ l/g 3/2 . The ratio of two predictions of the 
threshold velocity is Vj/vjj = g ^ 2 . Unfortunately, in 
our simulated systems segments sizes inside the OZ are 
of the order of goL,OR ~ 10- Thus, we cannot capture a 
significant difference between the threshold-velocity pre¬ 
dictions Vj and Vjj. In principle, brush systems with 
larger goL,OR values can be designed. However, for longer 
segment sizes, e.g. goL,OR ~ 100 chain entanglements 
will also come into play and introduce more complexities. 
Hence, throughout this work we will leave investigation 
of the threshold velocity v* for a separate work m , and 
refer v* only to distinguish between the linear (any seg¬ 
ment inside OZ is stretched less than its equilibrium size) 
and the non-linear force regimes (the segment end-to-end 
distance is more than its equilibrium size). 

As discussed above, at v > v* the segments inside the 
OZ are not Gaussian, and are stretched due to the rel¬ 
ative motion of brushes, as illustrated by the snapshots 
given in Fig. [7] At the velocity range v* < v < f max , 
where the maximum velocity that we considered here 
v max ~ 1 cr/r, both Sl and Sr decrease in a similar way 
with increasing velocity as shown in Fig. [4] Thus, Eq.[6]is 
still valid to explain the difference in the frictional forces. 

At around v max ~ 1 cr/r, regardless of the brush topol¬ 
ogy, the width of the OZs goes to unity 5l,r cr as 
shown Fig. [4j Hence, at v > r max , according to Eq. |6j 
the ratio of friction forces for linear and ring brushes is 
reduced to 


[l_ ~ 

Ir ^r 


(7) 


Indeed Fig.[6]shows that the ratio of binary-collision fac¬ 
tors for ring and linear Ql,r is even higher at high ve¬ 
locities and reaches VLr/VLr « 2. The limit of Sl,r cr 
for both ring and linear brushes at high velocities implies 
that there should be equal number of monomers inside 
the overlap volumes pm^L.R- Additionally, we know 
gR « 2 gL from Fig. |5| This is only possible if stretched 
segments of ring brushes form double-stranded confor¬ 
mations inside the OZ as can be seen in the snapshot 
given in Fig. [TJd. In the non-linear regime a segment 
of a ring brush with gR monomers inside OZ is highly 
stretched and adopts a double-stranded conformation. 
Due to double-stranding, each monomer of a segment in 
a ring brush has on average three neighbours from the 
same segment - two bonded and one nonbonded. This 
leads to smaller number of collisions with the monomers 
of the opposing brush. 

Another way of confirming that segments in the ring 
brushes are double-stranded inside the OZ is to calcu¬ 
late the fraction of “participating” chains, T l,r • Since 
linear chains has smaller number of monomers inside the 
OZ (gR < gR) larger number of chains should occupy the 
OZ to keep the monomer density uniform throughout 
the simulation box. Calculated values of T l,r according 


Ring brush 
@ Linear regime 



Linear brush; 

@ Linear regime 




FIG. 7: Snapshots of the individual chains from ring (a 
and b) and linear brush bilayer (c and d) simulations in 
the linear-force regimes (left column) and 
non-linear-force regimes (right columns). 
Polymerization degree of chains for linear brushes is 
N = 60 whereas for ring brushes is N = 120. Monomers 
of surrounding chains are rendered as points with colors 
referring to the top and bottom brushes. 


to previously described ” participating” chain description 
are shown in Fig. [8] with ratios given in the insets. For 
slow velocities (v < 10 _1 cr/r), which also corresponds 
to the linear-force regime, Tql is slightly larger than 
Tq r for all cases. On the other hand at high velocities 
v > v* the difference between linear and ring brushes 
increases and is consistent with the conditions gR > gL 
and 5l,r cr. Note that in Fig. [HJd, the increase in ratio 
is less pronounced. Possible reason of this is our criterion 
for chains participating in the OZ which underestimates 
looping segments since we check only the TVth monomer 
of the corresponding chain within the OZ. 

Finally, at ultra-high velocities v 1 cr/r, for which 
we cannot perform simulation since the thermostat can¬ 
not keep the system temperature uniform throughout the 
simulation box, one would expect that both Sl,r —>• cr 
and Ql,r 1- Since the grafted chains are completely 
inclined and gL,R 1, two opposing brushes can only in¬ 
teract via few monomers near the free edges of the grafted 
chains. This scenario indeed yields to a high-velocity lin¬ 
ear regime and /l/ /r —>• 1. 


CONCLUSION 

In this work we demonstrated using scaling arguments 
and MD simulations that friction forces of linear chains 
are higher than those of ring brushes for a broad range 
of velocities. For slow driving velocities (or shear rates), 
segments of a linear brush can penetrate through oppos¬ 
ing brush deeper compared to segments of ring brushes. 












9 



V [C/T] 

FIG. 8: The fraction of participating chains inside the 
overlap zone as a function of the plate velocity. Insets 
show the ratio of linear-brush data to ring-brush data. 

See text for details. The grafting densities are a) 

Pg = 0.25 a~ 2 and b) p^ = O.llcr -2 with p^ = 0.5p^. 


This is mainly due to the compactness of the segments 
in the ring brushes. The compactness also lessens the 
number of collisions between two opposing ring brushes 
and results in lower frictional forces between two rela¬ 
tively moving brushes. For large driving velocities (shear 
rates), segments inside the overlap zone between brushes 
are highly stretched regardless of the topology of chains. 
In the ring brushes, stretched segments form double- 
stranded conformations which reduces the number of col¬ 
lisions with the monomers of the opposing brush. 

We also observed that both friction forces and normal 
pressures in the linear chains are factor of two larger than 
those in ring brushes. This leads to equal kinetic friction 
coefficients p = f/p for both systems. However, the ef¬ 
fective viscosity rj e ff = f /v should be distinguishable for 
ring and linear chains. 

A more interesting situation can arise when entangled 
brushes are considered when chain segments can diffuse 
into the overlap volume with longer segments (around 
100 monomers per chain for the bead-spring model). 
Since the bulk melts of ring chains exhibit no relaxation 
plateau in their stress-relaxation moduli [26), 4fr, their 
frictional responses should be much lower than those of 
entangled linear brushes. We will consider this scenario 
in a future publication. 

We thank Edward J. Banigan and Ozan S. Sariyer for 
their careful readings of the manuscript. This research 


was supported by the Polish Ministry of Science and 
Higher Education (Iuventus Plus: IP2012 005072). 


[1] P G De Gennes. Polymers at an interface; a simpli¬ 
fied view. Advances in Colloid and Interface Science , 
27(3):189-209, 1987. 

[2] Jian Ping Gong. Friction and lubrication of hydrogels its 
richness and complexity. Soft Matter , 2(7):544, 2006. 

[3] B N J Persson. Sliding Friction: Physical Principles and 
Applications. Springer Berlin Heidelberg, 2000. 

[4] Neil Ayres. Polymer brushes: Applications in biomate¬ 
rials and nanotechnology. Polymer Chemistry , 1(6):769, 
2010. 

[5] A Halperin. Polymer Brushes that Resist Adsorption 
of Model Proteins: Design Parameters. Langmuir , 
15(7):2525-2533, March 1999. 

[6] Feng Zhou and Wilhelm T S Huck. Surface grafted poly¬ 
mer brushes as ideal building blocks for ?smart? surfaces. 
Physical Chemistry Chemical Physics , 8(33):3815, 2006. 

[7] Martien A Cohen Stuart, Wilhelm T S Huck, Jan Gen- 
zer, Marcus Muller, Christopher Ober, Manfred Stamm, 
Gleb B Sukhorukov, Igal Szleifer, Vladimir V Tsukruk, 
Marek Urban, Frangoise Winnik, Stefan Zauscher, Igor 
Luzinov, and Sergiy Minko. Emerging applications of 
stimuli-responsive polymer materials. Nature Materials , 
9(2): 101-113, November 2019. 

[8] Alexander Sidorenko, Sergiy Minko, Karin Schenk- 
Meuser, Heinz Duschner, and Manfred Stamm. Switch¬ 
ing of Polymer Brushes. Langmuir , 15(24):8349-8355, 
November 1999. 

[9] Philip Pincus. Colloid stabilization with grafted poly¬ 
electrolytes. Macromolecules , 24(10):2912-2919, 1991. 

[10] EPK Currie, W Norde, and M A Cohen Stuart. Tethered 
polymer chains: surface chemistry and their impact on 
colloidal and surface properties. Advances in Colloid and 
Interface Science , 100:205-265, 2003. 

[11] Uri Raviv, Suzanne Giasson, Nir Kampf, Jean-Frangois 
Gohy, Robert Jerome, and Jacob Klein. Lubrication by 
charged polymers. Nature , 425(6954): 163-165, Septem¬ 
ber 2003. 

[12] M Urbakh, J Klafter, D Gourdon, and J Israelachvili. 
The nonlinear nature of friction. Nature , 430(6999):525- 
528, 2004. 

[13] Jacob Klein, Eugenia Kumacheva, Diana Mahalu, Dvora 
Perahia, and Lewis J Fetters. Reduction of frictional 
forces between solid surfaces bearing polymer brushes. 
Nature , 370(6491):634-636, 1994. 

[14] J. Klein. Repair or Replacement A Joint Perspective. 
Science , 323(5910):47-48, 2009. 

[15] Meng Chen, Wuge H Briscoe, Steven P Armes, and Jacob 
Klein. Lubrication at physiological pressures by polyzwit- 
terionic brushes. Science , 323(5922): 1698-1701, March 
2009. 

[16] Bruno Zappone, Marina Ruths, George W Greene, Gre¬ 
gory D Jay, and Jacob N Israelachvili. Adsorption, Lu¬ 
brication, and Wear of Lubricin on Model Surfaces: Poly¬ 
mer Brush-Like Behavior of a Glycoprotein. Biophysical 
Journal , 92(5):1693-1708, March 2007. 

[17] Marcel Benz, Nianhuan Chen, and Jacob Israelachvili. 
Lubrication and wear properties of grafted polyelec- 






















10 


trolytes, hyaluronan and hylan, measured in the surface 
forces apparatus. Journal of Biomedical Materials Re¬ 
search , 71A(1):6-15, 2004. 

[18] N Maeda, N Chen, M Tirrell, and JN Israelachvili. Ad¬ 
hesion and friction mechanisms of polymer-on-polymer 
surfaces. Science , 297(5580):379, 2002. 

[19] Debby P Chang, Nehal I Abu-Lail, Farshid Guilak, Gre¬ 
gory D Jay, and Stefan Zauscher. Conformational Me¬ 
chanics, Adsorption, and Normal Force Interactions of 
Lubricin and Hyaluronic Acid on Model Surfaces . Lang¬ 
muir , 24(4): 1183-1193, February 2008. 

[20] Xavier Banquy, Joanna Burdynska, Dong Woog Lee, 
Krzysztof Matyjaszewski, and Jacob Israelachvili. Bioin¬ 
spired Bottle-Brush Polymer Exhibits Low Friction and 
Amontons-like Behavior. Journal of the American Chem¬ 
ical Society , 136(17):6199-6202, April 2014. 

[21] Debby P Chang, Nehal I Abu-Lail, Jeffrey M Coles, 
Farshid Guilak, Gregory D Jay, and Stefan Zauscher. 
Friction force microscopy of lubricin and hyaluronic acid 
between hydrophobic and hydrophilic surfaces. Soft Mat¬ 
ter , 5(18):3438, 2009. 

[22] Jeffrey M Coles, Debby P Chang, and Stefan Zauscher. 
Molecular mechanisms of aqueous boundary lubrication 
by mucinous glycoproteins. Current Opinion in Colloid 
& Interface Science , 15(6):406-416, December 2010. 

[23] J Dorier and A Stasiak. Modelling of crowded polymers 
elucidate effects of double-strand breaks in topological 
domains of bacterial chromosomes. Nucleic Acids Re¬ 
search , 41(14):6808-6815, August 2013. 

[24] Jonathan D Halverson, Won Bo Lee, Gary S Grest, 
Alexander Y Grosberg, and Kurt Kremer. Molecular dy¬ 
namics simulation study of nonconcatenated ring poly¬ 
mers in a melt. I. Statics. The Journal of Chemical 
Physics , 134(20) :204904, 2011. 

[25] M E Cates and J M Deutsch. Conjectures on the statistics 
of ring polymers. Journal de Physique , 47(12) :2121, 1986. 

[26] Jonathan D Halverson, Won Bo Lee, Gary S Grest, 
Alexander Y Grosberg, and Kurt Kremer. Molecular dy¬ 
namics simulation study of nonconcatenated ring poly¬ 
mers in a melt. II. Dynamics. The Journal of Chemical 
Physics , 134(20):204905, 2011. 

[27] Yangpeng Ou, Jeffrey Sokoloff, and Mark Stevens. Com¬ 
parison of the kinetic friction of planar neutral and poly¬ 
electrolyte polymer brushes using molecular dynamics 
simulations. Physical Review E, 85(1), January 2012. 

[28] L Spirin and T Kreer. Strongly Compressed Poly¬ 
electrolyte Brushes under Shear. ACS Macro Letters , 
2(l):63-66, January 2013. 

[29] A Galuschko, L Spirin, T Kreer, A Johner, C Pastorino, 
J Wittmer, and J Baschnagel. Frictional Forces between 
Strongly Compressed, Nonentangled Polymer Brushes: 
Molecular Dynamics Simulations and Scaling Theory. 
Langmuir , 26(9):6418-6429, May 2010. 

[30] Jan-Michael Y Carrillo, Daniel Russano, and Andrey V 
Dobrynin. Friction between Brush Layers of Charged 
and Neutral Bottle-Brush Macromolecules. Molecular 
Dynamics Simulations. Langmuir , 27(23): 14599-14608, 


December 2011. 

[31] Gary S Grest and Michael Murat. Structure of grafted 
polymeric brushes in solvents of varying quality: a molec¬ 
ular dynamics study. Macromolecules , 26(12):3108-3117, 
1993. 

[32] Benoit Liberelle and Suzanne Giasson. Friction and 
Normal Interaction Forces between Irreversibly At¬ 
tached Weakly Charged Polymer Brushes . Langmuir , 
24(4):1550-1559, February 2008. 

[33] K Binder and A Milchev. Polymer brushes on flat and 
curved surfaces: How computer simulations can help to 
test theories and to interpret experiments. Journal of 
Polymer Science Part B: Polymer Physics , 50(22): 1515- 
1555, September 2012. 

[34] TA Witten, L Leibler, and PA Pincus. Stress relaxation 
in the lamellar copolymer mesophase. Macromolecules , 
23(3):824-829, 1990. 

[35] Alexander N Semenov. Rheology of Polymer Brushes: 
Rouse Model. Langmuir , ll(9):3560-3564, September 
1995. 

[36] J F Joanny. Lubrication by molten polymer brushes. 
Langmuir , 8(3):989-995, 1992. 

[37] Jacob Klein. Shear, friction, and lubrication forces be¬ 
tween polymer-bearing surfaces. Annual review of mate¬ 
rials science , 26(1):581-612, 1996. 

[38] D Reith, A Milchev, P Virnau, and K Binder. Anoma¬ 
lous structure and scaling of ring polymer brushes. EPL 
(Europhysics Letters ), 95(2):28003, June 2011. 

[39] Fang Yin, Dmitry Bedrov, Grant D Smith, and S Michael 
Kilbey. A Langevin dynamics simulation study of the tri¬ 
bology of polymer loop brushes. The Journal of Chemical 
Physics , 127(8):084910, 2007. 

[40] Kurt Kremer and Gary S Grest. Dynamics of entangled 
linear polymer melts: A molecular-dynamics simulation. 
The Journal of Chemical Physics , 92(8):5057, 1990. 

[41] C Pastorino, T Kreer, M Muller, and K Binder. Compar¬ 
ison of dissipative particle dynamics and Langevin ther¬ 
mostats for out-of-equilibrium simulations of polymeric 
systems. Physical Review E, 76(2):026706, August 2007. 

[42] Steve Plimpton. Fast parallel algorithms for short-range 
molecular dynamics. Journal of Computational Physics , 
117(1): 1—19, 1995. 

[43] W. Humphrey, A. Dalke, and K. Schulten. VMD pro¬ 
gram. J. Molec. Graphics , 14(1):33, 1996. 

[44] Florent Goujon, Aziz Ghoufi, Patrice Malfreyt, and Do¬ 
minic J Tildesley. Frictional forces in polyelectrolyte 
brushes: effects of sliding velocity, solvent quality and 
salt. Soft Matter , 8(17):4635-4644, 2012. 

[45] A. Erbas, Zhulina E., and M. Rubinstein. Viscous friction 
of linear polymer brushes, in preparation. 

[46] M Kapnistos, M Lang, D Vlassopoulos, W Pyckhout- 
Hintzen, D Richter, D Cho, T Chang, and M Rubinstein. 
Unexpected power-law stress relaxation of entangled ring 
polymers. Nature Materials , 7(12):997-1002, December 
2008. 



